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S e presentan dos metodos que se pueden 
utilizar como alternativas al GWAS para 
localizar QTLs en ganado. El primero se 
centra en el analisis de secuencias de homo- 
cigosis (runs of homozygosity, ROH) en gru- 
pos de animales (p.ej. cases y controles); el 
segundo es un metodo de remuestreo basado 
en la frecuencia de inclusion de los SNP en 
modelos predictivos. Las ROH se aplicaron a 
la identificacion de regiones del genoma aso- 
ciadas a trastornos reproduct ivos en vacas de 
raza Prisona, y el metodo de remuestreo a la 
deteccion de portadores del haplotipo BH2 en 
vacas de raza Pardo Alpina. Estos metodos 
alternatives pueden complementar el GWAS 
estandar en la localizacion de QTLs para car- 
acteres de interes en el genoma de animales 
domesticos. 

[English abstract] 

eyond GWAS: alternatives to localize 
QTLs in farm animals. Two methods 
that could be used for QTL mapping 
as alternatives to standard GWAS are pre¬ 
sented. The first relies on the differential 
frequency of runs of homozygosity (ROH) in 
groups of animals (e.g. cases and controls), 
while the second stems from resampling tech¬ 
niques used for the prediction of carriers of 
a mutation, and is based on the frequency 


of inclusion of polymorphisms (SNP) in the 
predictive model. ROH were applied to the 
detection of reproductive diseases in Holstein- 
Priesian cattle, while resampling was applied 
to the detection of carriers of the BH2 haplo- 
type in Brown Swiss cattle. These alternative 
approaches may complement GWAS analyses 
in localizing more accurately QTLs for traits 
of interest in livestock. 

Introduccion 

El descubrimiento de regiones del genoma asociadas 
a caracteres de interes zootecnico (p.ej. resistencia 
a enfermedades) es una importante etapa prelimi- 
nar para posibles aplicaciones de la genomica a la 
cria de animales domesticos como la identificacion de 
portadores de mutaciones favorables/desfavorables o 
la seleccion asistida por marcadores. Esto es ahora 
posible gracias a la gran cantidad de datos proce- 
dentes de las nuevas tecnicas de secuenciacion masiva 
(NGS: next generation sequencing) en forma p.ej. de 
genotipados de alta densidad de SNPs y de secuen¬ 
cias de genoma completos. Los estudios de asociacion 
de genoma completo (GWAS) se han convertido 
por lo tanto en una tecnica estandar para escanear 
el genoma en busca de polimorfismos asociados al 
fenotipo analizado. Sin embargo, los metodos de 
GWAS presentan algunas limitaciones: analizan gen- 
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eralmente un SNP a la vez, son susceptibles de pro- 
ducir resultados espurios y, salvo en caso de senales 
de asociacion muy patentes, pueden resultar de difcil 
interpretacion (McCarthy et al. (2008)) . En este 


arti'culo se presentan dos alternativas al GWAS para 
localizar QTLs en ganado: el analisis de secuencias 
de homocigosis (“runs of homozygosity”, ROH) y 
un metodo de remuestreo basado en la frecuencia de 
inclusion de los SNP en modelos predictivos. 


Material y Metodos 

Para ilustrar las diferentes propiedades de los 
metodos basados en GWAS, ROH y remuestreo para 
localizar QTLs, se utilizaron los resultados de tres 
estudios previos con bovinos lecheros de raza Frisona 
y Pardo Alpina. En particular: el GWAS se baso en 
17 vacas afectadas por el desplazamiento abomasal 


izquierdo y 47 controles de raza Frisona (Biffani et al 


(2014)); el analisis de ROH en 163 vacas afectadas 
por trastornos reproductivos y 295 controles de raza 


Frisona (Biscarini et al. (2014b); y el metodo de re¬ 


muestreo en 513 portadores del haplotipo BH2 en 
BTA19 y 3132 no portadores de raza Pardo Alpina 


(Biffani et al. (2015)). Todos los animales se genoti- 
paron con el chip bovino 50k de 54 001 SNPs. 


ficaron entonces las ROH presentes en la poblacion 
de vacas bajo estudio a traves de la estimacion de la 
homocigosis de los SNP en “ventanas deslizantes” de 
1000 kb. Se permitieron como maximo 5 genotipos 
faltantes y 1 genotipo heterocigoto para considerar 
la secuencia como ROH. La frecuencia relativa de las 
ROH en casos y controles indica regiones del genoma 
asociadas con el fenotipo. 

Remuestreo de SNPs en modelos predictivos 

Para la identificacion de portadores del haplotipo 
BH2 en vacas de raza Pardo Alpina se aplico un 
procedimiento basado en la reduccion progresiva del 
mimero de SNP incluidos en el modelo y en el analisis 
discriminante lineal (LDA). Para cada umbral de 
SNP utilizados (2.5%, 10%, 15%, 30%, 50% y 100% 
sobre el total de SNP) se seleccionaron aquellos mas 
predictivos a traves del metodo de Best Subset Se¬ 
lection (BSS) y se utilizaron los SNP seleccionados 
para clasificar los individuos en portadores o no del 
haplotipo. El procedimiento se repitio 1000 veces 
(validacion cruzada de 10 pliegues, 100 replicas) por 
cada umbral, y se obtuvo la grafica de la frecuencia 
de inclusion de los SNP en el modelo predictivo en 
funcion de su posicion en BTA19. De esta manera 
se pudo localizar la mutacion en el genoma. 


GWAS : El estudio de asociacion pangenomico 
(GWAS) para el desplazamiento abomasal izquierdo 
se llevo a cabo siguiendo el procedimiento 
GRAMMAR-GC descrito por Aulchenko et ah 
(2007): primero, se ajusto un modelo de re- 
gresion logfstica para datos binarios (caso-control) 
incluyendo los efectos sistematicos de granja, orden 
de parto y el efecto poligenico aleatorio, y posterior- 
mente se utilizaron los residuos de este modelo en 
una regresion lineal para estimar el efecto de cada 
SNP. 


Analisis de secuencias de homocigosis (ROH) 

: Bajo la hipotesis de que las enfermedades comple- 
jas presentan una componente gentica constituida 
por muchas variantes recesivas distribuidas a lo largo 
del genoma, se analizaron las secuencias de homo¬ 
cigosis (ROH) para identificar regiones asociadas a 
trastornos reproductivos en vacas lecheras. Las ROH 
se definen como secuencias de genotipos homocigotos 
adyacentes, que reflejan la transmision de haplotipos 
identicos por parte de progenitores comunes. En 
lugar de centrarse en un locus, las ROH consideran 
tambien al entorno, teniendo asf en cuenta el cambio 
de frecuencia de SNP cercanos que esten en desequi- 
librio de ligamiento con el locus analizado. Se identi- 


Resultados y Discusion 

En el estudio de asociacion se detectaron SNPs aso- 
ciados con el desplazamiento abomasal izquierdo en 
BTA12. La asociacion mas significativa se halla a 
una distancia de tan solo 20 kb del gen SLITRK5. 
Este gen pertenece a la familia genica SLITRK, impli- 
cada en procesos neurologicos que ya se demostraron 
asociados a la patogenesis del desplazamiento abo¬ 
masal en vacuno. La Figura[^ present a el Manhattan 
plot del analisis GWAS. Se puede apreciar la falta 
de senales claras de asociacion y la aparicion de nu- 
merosas asociaciones, potencialmente espurias, que 
resultan mas o menos “significativas” y distribuidas 
a lo largo del genoma. 

La interpretacion de los resultados puede ser arbi- 
traria, lo que ilustra algunas de las limitaciones del 
GWAS. El analisis de las secuencias de homocigosis 
(ROH) permite ampliar la mirada al entorno de cada 
SNP, llevando a interpretaciones mas robustas de los 
resultados: la aparicion de un haplotipo homocigoto 
en casos pero no en controles puede ser mas infor- 
mativa que un SNP asociado individualmente. La 
Figuraj^ ilustra el caso de trastornos reproductivos 
en vacas lecheras en BTA15. El analisis de ROH es 
un metodo falto de un modelo estadfstico explfcito. 
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Figure 1: Manhattan plot de los resultados (log{pjvalue)) del estudio de asociacion entre SNP y 

desplazamiento del abomaso en vacas de raza Frisona Holandesa. La Unea roja indica el umbral de 

significacion: —/o5(0.0001) = 4 


y esto puede causar problemas a la hora, p.ej., de 
probar la significancia de la asociacion o de incluir 
efectos sistematicos en el analisis. 


Sin embargo, hay maneras de resolver estas difi- 
cultades. La significancia estadfstica se puede testar 
comparando la homocigosis media (o la frecuencia 
de las ROH) en casos y controles con un test-t, o 
aplicando el concepto clinico de “no inferioridad” 

. Se realizan dos analisis, 
ROH y GWAS, y se calcula el “false discovery rate” 
(FDR) para los dos: la hipotesis nula es que el FDR 
es mayor en ROH que en GWAS (ROH inferior 
al GWAS), la hipotesis alternativa es que los dos 
metodos son equivalentes (mismo FDR medio). Los 
efectos sistematicos se pueden incluir utilizando un 
modelo estadistico para analizar las ROH 
; otra alternativa es estratificar el 
de ROH segiin las clases del efecto sistematico (ej. 
machos/hembras), o aplicarlo a los residuos de un 
modelo lineal. Mas detalles se encuentran en IBis-l 


(2012 


(Pollott 
analisis 


(D’Agostino et al. (2003 


carini et al. (2014a|b ). 


Otra opcion para localizar QTLs es utilizar los 
resultados de tecnicas de remuestreo. La Figuraj^ 
presenta la frecuencia de inclusion de los SNP en el 
modelo predictivo para portadores del haplotipo BH2 
en funcion de la posicion en BTA19. Las tecnicas de 
remuestreo permiten analizar la variabilidad alrede- 
dor del parametro estimado y se caracterizan por 
tener mas poder estadistico y localizar QTLs de man- 
era mas robusta e independiente de los valores p y 
del ensayo de frecuencia a cada locus. Los metodos 


propuestos pueden ser utilizados para complemen- 
tar el GWAS estandar en la localizacion de QTLs 
para caracteres de interes en el genoma de animales 
domesticos. 
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Figure 2: Valores de homocigosis a lo largo del BTA 15 en casos (rojo) y controles (azul) para enfermedades 

reproductivas. 
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Figure 3: Frecuencia de inclusion de los SNP en el modelo predictivo (en cada de las 1000 replicas por 

umbral) en funcion de la posicion en BTA 19. 
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